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Abstract 

The Loewner equation describes the time development of an ana- 
lytic map into the upper half of the complex plane in the presence of 
a "forcing" , a defined singularity moving around the real axis. The 
applications of this equation use the trace, the locus of singularities in 
the upper half plane. This note discusses the structure of the trace for 
the case in which the forcing function, £(£), is proportional to (— 1)@ 
with (3 in the interval (0,1/2). In this case, the trace is a simple 
curve, 7(t), which touches the real axis twice. It is computed by using 
matched asymptotic analysis to compute the trajectory of the Loewner 
evolution in the neighborhood of the singularity, and then assuming a 
smooth mapping of these trajectories away from the singularity. Near 
the t — singularity, the trace has a shape given by 

(Re ( 7 (t) - 7 (0)) 1 - /3 « (Aim 7 (t)) = 0(£(*)) 1 " /3 

A numerical calculation of the trace provides support for the asymp- 
totic theory. 
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1 Introduction 



1.1 Defining Loewner evolution 

0. Schramm has pioneered a new mathematical approach to the geometry 
of fractal and/or critical objects in two dimension. Many different applica- 
tions have been found to such problems as percolation, random walks, and 
critical phenomena 1 . The approach has been called Schramm-Loewner evo- 
lution (SLE) because it is based upon the earlier work of C. Loewner [2] on 
the relation between two-dimensional geometry and differential equations 
in the complex plane. This approaches is based upon the relatively simple 
equation 

dgt __ 2 

dt 9t - m uaj 

containing the forcing function, a real function of a real variable. Equa- 
tion^] is solved with the initial condition 

9to = z (lb) 

where z is a complex variable that lies in the upper half of the complex 
plane, thereby generating the function gt(z). For each value of t > to, this 
function can be thought of as a mapping of the form w = gt(z) which takes 
some connected subset of the upper-half z-plane, 1Z, into the entire region 
above the real axis of the w-plane, U. It also maps the curve bounding 1Z, 
called r, onto the real line of the u>-plane. 

The Loewner equation continually generates new singularities of gt(z). 
At each time, s, a new singularity is generated at the point, z = j(s), such 
that 

9s(l(s))=m- (2) 

As time goes on, points in the upper half-plane are continually added to 2 . 
7 and T. The set T is a geometrical object which has evoked considerable 
interest, particularly in the case in which the forcing, is a Brownian 
process (Zj. 

The function 7(i) has some very remarkable geometrical properties, all 
derived from the nature of the singularities in £(t). For the purposes of this 
paper, the most important such property is the non-self-touching behavior 
of r. If has all its singularities weaker than (t c — t) 1 / 2 (for example 

1 A bibliography can be found in [7j. 

2 If 7 is a simple curve that intersects the real axis once, Y is the union of that curve 
and the real axis. 
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if it has a derivative everywhere) then 7(i) is a simple curve which never 
touches itself. Conversely if has a singularity stronger than the one half 
power, e.g. (t c — t) 13 with < (3 < 1/2 then T will intersect or touch itself 
[S] . These touchings (or non-touchings) are crucial to the descriptions of the 
critical processes which are the main applications of SLE. One example of 
such a touching trace is shown in figure ^ It rises on the left as a square 
root of time, but is there non-singular in its shape. On the right is a singular 
osculation with the real axis. This paper is devoted to finding the form of 
this singularity. 

1.2 Orbits 

To calculate the trace we need to find the orbits or trajectories produced by 
equation Hal Thus we look at gt{z) while varying t with z held fixed. One 
family of such orbits runs from the trace to the real axis. It is formed by 
setting z = j(s) where s is a real variable with the meaning of a time. We 
demand s > to- Then the orbits we consider are the curves traced out by 
9t(l( s )) when t is varied through the interval [ioi 5 ]- At the one end, the 
curve hits points on the trace, since 

9t (l(s))=l(s), (3) 

On the other end, the curve intersects the real axis (see equation |5J) • The 
intersection is orthogonal if is non-singular at t = s. To find the trace 
[3] use equation|2]to provide initial conditions for all the trajectories gti'y(s)) 
and run the evolution, equation > backward from time s to time to- For 
each s, the solution at that time determines j(s) via equation^ One family 
of such orbits is depicted in figure El 

This kind of analysis is familiar from the theory of the motion of "passive 
scalars" (Hj. Take equation Hal to describe the motion of a point in the 
w— plane driven by a velocity field 2/{w — £). At first sight, this appears to 
be just the motion engendered by a point source or sink at w = £. However 
that interpretation is wrong. A point sink at £ would generate the field 
—2/(w — £) where the bar stands for complex conjugation. Figure [21 shows 
the difference between the two flow patterns. The conjugated field has an 
outflow per unit time of fluid occupying one unit of area. Thus after a 
time, an entire area of the fluid has escaped. The actual flow of equation Hal 
directly allows a much more restricted escape. If the forcing is sufficiently 
smooth, everywhere smoother than Hoelder index 1/2, the escaped points 
form a simple curve which intersects the real axis at £(io)- 
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Figure 1: This figure shows the trace for the case of power law forcing for 
(3 = 1/4, k = 8 and irj = —1,0. The right hand side of the trace curve is 
a nonsingular orthogonal intersection with the real axis. On the left the 
singularity in the forcing has generated a singular osculation with the axis. 

1.3 Outline of paper 

This note describes the singular properties of equation when the forcing 
is singular and has the form 

m = 2^{-tf (4) 

for the case in which (3 is in the interval (0,1/2) and to is smaller than 
zero. To understand the various different t = singularities generated by 
the Loewner evolution, we concentrate upon the trace, 7(i). All other sin- 
gularities are derived from the ones in j(t). 

Both trajectories and trace can be studied analytically [S] for (3 = 0, 1/2, 
and 1. In this paper we study singularities for (3 in the interval (0, 1/2) by 
performing an asymptotic analysis of the properties of trajectories, gt(z), 
for t and z in the neighborhood of the singularity. 

The next section of this note is devoted to describing the exact solutions. 



4 




\ \ 

c d 

Figure 2: These four panels show a "fluid" in grey. The grey arrow points 
to a singularity or leak on the boundary. In panel a, the flow is the one 
produced by a point sink; in b, the flow pattern is the one used in this paper 
(see equation ED) in which there is downflow in the vertical and outflow in 
the horizontal direction. The white regions in the next two panels show the 
regions containing the fluid which will leave the system through the leak. 

and the orbits which arise from them. Results for (3 = 1 are particularly 
needed because they are used in the subsequent small-/? analysis. Section |21 
turns to the situation for < (3 < 1/2, which is studied via matched asymp- 
totic analysis 6 . This part looks at the behavior of the critical trajectory, 
the one which passes through z = 0, and at the linear stability analysis 
about that trajectory. The critical trajectory has the unique property that 
it intersects the real axis twice. Section 13.21 is devoted to describing the 
small z— trajectories near their intersection with the real axis. The final 
section connects the behavior of these trajectories to the linear analysis of 
section 3 and thereby finds the actual singularity in the trace. The trace is 
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found asymptotically when — to is a positive number much smaller than one. 
Away from the singularity at time zero, the time development is described 
by smooth mappings, and consequently we argue that leading order shape 
of the trace is universal and applies, for example, to larger values of — to- 

2 Exact Analysis 

2.1 The case (3 = 0. 

Pick £(t) = 0. Then the solution for g is simply 

g t (z) = z 2 + 4(t - t ) (5) 

and the trace lies on the imaginary axis, having j(s) = 2i^J s — to. As a 
result the orbit leading up to the singularity is 

9t(l(s)) = 2iV~s - t for t < s < t (6) 

In this case, for all values of s the singular trajectories gt(j(s)) follow the 
imaginary axis. 

2.2 The case (3 = 1. 

Next, consider the case of a linear forcing: 

m = to - 1 (7) 

which was fully solved in reference [5]. The linear solution 3 depends upon 
the function F which is defined by 

F(z) = -z + 2\n(l + z/2) (8) 

This paper will make use of the functional inverse of F, i.e. the function 
with complex x in the upper half plane, which obeys F~ 1 (F(z)) = z. 
For x with small magnitudes, 

F^ix) = 2iy/x~ - 2x/3 + —x 3 / 2 + • • • (9a) 
18 

3 Here we use a convention for the sign of t in equation |7] different from the one used 

in ref 0. This results in sign changes in equation |H| compared with the corresponding 

equation in that paper. 
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The branch is picked so that the square root is positive when x approaches 
the real axis. If the real part of x is large then the corresponding branch 
gives 

F~ l (x) = -x + 2pn(s) + m-hx2}[l - 2/x] - 4/x + • • ■ (9b) 

The solution for g is given in terms of this known function, F, by 

F(g t (z)-m=F(z)+m (10) 

The trace is given by setting z = j(t) with 7 being picked so that gt(l(t)) = 
£(i). In this situation, this identification makes the left hand side of equation 
ITUl equal to -F(O) = 0, and consequently the trace obeys 

F( 7 (s)) = -as) = s-t 

with the final consequence 

7 ( fl ) = *-*(-£(,)). (ii) 

On the other hand, the critical trajectory is obtained by setting z = 'j(s) in 
equation El which then reads 

F(g t (i(s)) - m) = ni(s)) + m = s-t 

This equation may then be solved to find the orbit in terms of the trace. 

gt(7(s))=m + F~ l (s-t) (12) 

Figure shows the behavior of the trace and some of the trajectories. The 
trajectories connect the trace-points with the corresponding singular points 
on the real axis. In fact, of course, the trace is calculated as the locus of end- 
points of the trajectories. At the intersections of the orbits with the real axis, 
the curves meet the axis orthogonally. (Also see equation I13bl below.) This 
behavior is required whenever £(t) is non-singular at the intersection. The 
points plotted on the orbits are all equally spaced in t. The large spacings 
on the figure near the intersections reflect the square-root singularities at 
these points. The figure hints that even for large s values the trajectories 
remain bounded below y = 2ir, and only grow slowly in their overall length. 
The exact solution (see equation I14bl below) shows that this growth is, in 
fact, logarithmic in s. 
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Figure 3: This figure includes the trace and three trajectories for the case 
of linear forcing. Specifically, we take £(t) = (5 — t) and to = 0- The solid 
line at the top is the trace. The symbols define trajectories with, counting 
from the left, s = 5, 20, and 40. 



The results in figure |3] can be seen more precisely by calculating the 
asymptotic properties of trace and orbit from the exact solution. For s just 
bigger than to , equation implies that the trace has the value 

7(a) = 2iy/s=to~-2(s-t ) /% + ■■■ (13a) 

while from equation^] as t approaches s from below (and s > to), i.e. just 
before the trace hits the singularity, the orbit obeys 

gt(l(s)) = t - t + 2iV^t - 2(s - t)/3 + • • • (13b) 

In the opposite limit, i.e. as time, s, goes to infinity, the trace obeys 

7 ( s ) = s + t + 2 ln[(s - t )/2] + 2ni + 0( ln(s ~/ o) ) (14a) 

s - 1 

while, for large s — t the orbit approaches the limit 

9t{l{s)) =to-s + 2[ln[(s - t)/2] + 2m][l - 2/{t -*)] + ••• (14b) 
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Eauation ll4al shows as s becomes very large, the trace moves with almost the 
same velocity as the point of forcing, £(s), but it always lags slightly behind 
it. Asymptotically, the separation between trace and £ goes to a constant, 
2tt, in the imaginary direction and grows only logarithmically in the real 
direction. Correspondingly, equation I14bl shows that the orbit covers this 
separation in the time between to and s. 

2.3 The case (3 = 1/2. 

So far we have seen no intersections or touching of trajectories and the 
real axis, except for the trivial take-off at to- With a forcing of the form 
lyj k(1 — t), and k > 4, there is an exact solution 7 which exhibits a non- 
trivial intersection with the real axis. The intersection is approached by a 
trace in the form of straight lines at an angle to the real axis. The angle is 
expressed in terms of the parameters [7], y± 

V± = ± - K (15) 
The angle is given as eft where 

4> = ny-/y+. (16) 

This trace and its comparison with the exact solution are shown in figure 
El for the case in which k = 16. At s = 0, the value of 7(5) is equal to £(s) 
and is hence 8. This value is the point on the far right of the trace curve. 
From this point, as s increases, the imaginary part of the trace grows and 
then falls away. Over the entire range of the plot, the real part of the trace 
is roughly proportional to £(s). Finally as s approaches the singularity at 1 
the imaginary part of the trace decreases and goes to zero proportionally to 
£(s), In this limit, the real part also varies as £. Thus as s approaches 1 the 
trace has the behavior 

1 {s)-y.=Ce i H{s). (17) 

where <f> is given by equation El C is real and positive, and y_ is the value 
of the trace at zero found from the exact solution. 

Figure H] was calculated numerically, rather than from using the exact 
solution. Since there is a singularity at s = 1, the calculation could not be 
carried all the way down to s = 1. However the calculation did continue 
quite close to the exact intersection point at 7(1) = y_ 0.53, specifi- 
cally to a time within 10 -5 of the critical time. The calculation is rather 
straightforward and is described in the appendix. Since the exact solution is 
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Figure 4: This figure compares the trace for a square root singularity and 
k = 16 with its small-s fit. The numerical data is the points; the line is the 
fit with slope and intercept chosen to be the right theoretical values. The 
fit looks and is excellent. 

available, a comparison to the exact solution can be made, giving a relative 
accuracy of 3 parts in 10 s at this closest trace point to the singularity, and 
better accuracy further away. This curve also shows that the numerically 
calculated result is fit over a quite reasonable range by the near-singular 
asymptotics. 

3 The Critical Trajectory 

Here we discuss the singularities in Loewner evolution for the case in which 
the forcing has a singularity with index (5 at the time t = 0. As we shall 
see, in the studied interval, < (3 < 1/2, for values of t near zero the 
critical trajectory has gt much smaller than Hence this trajectory and 
nearby ones can be studied by expanding in the g on the right hand side 
of equation Hal This expansion produces an analytically tractable problem 
which is studied in this section. 
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3.1 Behavior of singular trajectory 

We solve equation Hal going backward in time from the "initial" condition 

3o(7(0)) =0 (18) 

We notice that £ is changing quite rapidly in the neighborhood of t = 
so that we might expect that ^(7(0)) is much smaller than £(i) in this 
neighborhood. Thus we might try an expansion of equation Hal of the form 

^(7(0)) = -2/f - 2 5i ( 7 (0))/£ 2 - 2(g t ( 1 (0))) 2 /e + ■■■ (19) 
This equation can be solved iteratively with the lowest order solution being 

ft(7(0))= +••• 

and the next order 

*M0)) = ^L + _i_l_ + ... po. 

This trajectory is the one which causes the trace to hit the real axis. It has 
this effect because it remains real. This reality is the expected behavior for 
a trajectory which describes a point at which the trace touches or intersects 
itself. On the other hand, all the other trajectories pop off the real axis as 
soon as —t increases beyond its initial value, —s. 

Convergence of this series requires (3 to be less than one half so that 
exponents of the form n — (2n + \)(3 do not change sign. 

Note also that our initial assumption that gt <C is quite true for this 
trajectory since 

{-tf-P « {-tf 

for small values of —t and (3 in the interval (0, 1/2). 

Thus equation l20l is accurate when —t <C 1. If it is also true that —to *C 1. 
we can find find the g-value at the initial time to < and hence find the 
trace at s = as 

7(0) = 9t (7(0)) = — p + {1 _ m _ 3p) (21) 

This then tell us the position of the end of the trace, but nothing about the 
shape of the curve 7(s). Further, we do not know the trace even at s = 
when —to is not small. 
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3.2 Near-critical behavior: Linear deviations from critical 
trajectory 

To see some of the properties of a trajectory near the critical one, take s to 
be small and write 

9t(l(s)) = 9t(l(0)) + A t . (22) 

Then the deviation At must in some sense be small. In a linear analysis, it 
obeys 

lAt = -2At/[g t ( 7 (0))-m? (23) 
with the initial condition that 

A s = 0.(700) - 5,(7(0)) » (sf - 2( ~y (24) 

However, we cannot use this initial condition here. Because the first term 
dominates the right hand side of equation 1241 At is not a small perturbation 
near t = s and we cannot employ equation 231 in the region. Thus we cannot 
make direct use of the initial data of equation [23] in our analysis of equation 

[23 

On the other hand, for larger values of —t, i.e. specifically for some 
negative t in the range. — s <C — t <C 1 equation [221 does apply. Under the 
stated conditions, g is much smaller than £ so that the equation takes the 
simple form 

|lnA t = -2[£(t)]- 2 
which then has the solution 

9t(l(s)) = 5,(7(0)) + A r exp[2i T ^- - 2 ^^ ] 

where T is some reference time in the interval — s <C — T <C 1. Since T has 
no physical meaning, one can absorb the factor involving T into At and get 
a result of the form 

g t ( 7 (s))=gt(7(0)) + Aexp[2 Yr^ ] (25) 

Equation [23 is accurate only when one can expand in the correction 
term in A. If — is bigger than or of the order of one and -t< 1, the 
explicit i-dependence in the exponent of Equation [23 does not matter, and 
the equation is accurate whenever 

|A| < \((t)\ and -i>-s (26) 
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To complete this analysis and find our trace, we have to extrapolate to 
values of t closer to s, and then used the boundary condition of equation 
1241 This analysis is the job of the next section. The final section brings the 
various results together. 

4 Near-critical trajectories close to their intersec- 
tion with the real axis 

In the previous section, we looked at situations in which £(t) 3> gt- Now 
we wish to focus upon the orbits' intersections with the real axis. At inter- 
section, = gt; near intersection they are close to one another. So an 
expansion of the right hand side of equation Hal in gt will not work. 

On the other hand, near enough to the intersection we can expand £(t) 
in a power series around t = s. Because the linear case is exactly solvable, 
we can expect to get accurate solutions (albeit with only a limited range of 
validity) by using this approach. 

The trace, 7(5), has been defined so that the trajectory g t ('f(s)) intersects 
the real axis at t = s. Take s to be a small negative number. In the 
neighborhood of the intersection the forcing has the value 

£(s-t) « £ l (s -t) = <5 + tA for < r < -s (27) 

where the constants 5 and A are respectively very small, 

5 = 2^{-sf. (28a) 

and very large 

A = 2^f3{-sf- 1 (28b) 

4.1 Solution of equation with linearized forcing. 

Equation EH gives the solution when the forcing is just —t. Our case is 
slightly different for two reasons. First the forcing is not just —t but a linear 
function of t. Secondly because the solution derived in this section only 
applies for small r, we cannot directly apply a boundary condition in z. To 
make these changes, we use the invariances in ref (Hj to see that the solution 
takes the form: 

F{A[g t (z)-( L (t)\) =c(z)+A£ L (t) 

where F is given by equation |HJ c{z) is a "constant" of integration, and 5 
and A are both independent of r. To set the constant of integration for the 
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trajectory, note that at the intersection of the trajectory with the real axis, 
£ = 5, while the argument of F vanishes and so does F itself . Thus we find 

F{A[g^ T { 1 {s))-5-Ar]) = A 2 r 

which then has the solution 

g s - T ] ( 7 («)) = 5 + At + A^F^^t) (29) 

This result can be stated explicitly in two limits. From equation we find 

g s - T (l(s)) =S + Ar/3 + 2i^ for < A 2 r < 1 (30a) 

The opposite limit comes from equation I9bl which then gives us the limiting 
form 

9s-t(j(s)) = 5 + A -1 [27rf + ln(A 2 r/2)] for A" 2 < r < -s (30b) 

Equation I30bl is the main result of this section. This equation enables 
us to see that for t much larger than |s|, to leading order, g contains 

• a large real s— dependent term (5 = 2 v / k(— s) 13 ) 

• a much smaller imaginary s— dependent term, and 

• a real term with a weak t— dependence. 

The last dependence is weak both because the coefficient is small and because 
a logarithmic dependence is a weak one. 

5 Connections 
5.1 A small difficulty 

This paper has concerned itself with two different expansions, the expan- 
sion about the critical trajectory of section |3] and the linear expansion of 
£ in section 21 In the language of matched asymptotics |S] , the latter is 
a sort of inner expansion, the former a sort of outer one. So in the usual 
approach, we determine the overlap region of the two expansions, determine 
the parameters on one in terms of those on the other side, and declare the 
problem solved. The linear expansion for gt(j(s)) is acceptable so long as 
t = t — s <S s. Conversely the expansion about the critical trajectory is 
sensible when (see equation ) the deviation parameter A from the crit- 
ical trajectory obeys |A| <C £(i). The most reasonable guess for A, which 



14 



is a deviation from criticality for 54(7(5)), it to choose it to be the biggest 
constant term in g as it was determined in the linear expansion, namely 
& = £( s ) ~ ( — ■ Now that result presents a problem. Our deviation from 
criticality argument is correct only if r >> —s, exactly the opposite of the 
range permitted by the expansion argument. We now conclude that there 
is no overlap of regions of validity of the two expansions, and consequently 
we must do better to match up our asymptotics. 

5.2 A better expansion 

This section will make heavy use of the result of equation I30bl which says 
that, in an appropriate range, 54(7(5)) is almost independent of t. The 
starting point will be a consequence of equation Hal that 



Take T to lie in the asymptotic large-T range of the linear expansion (see 
equation I30b|) and substitute the known form of g in this range for gx- Thus 
our integral expression for g becomes 



5t(7(s)) = 5 + A- 1 [ln(A 2 (s - T)/2) + 2vri] - f T 2 du/[£(u) - 9u(l(s))] 



To understand the integral equation, notice that the integral has a loga- 
rithmic singularity near its lower cutoff at T. This singularity arises because 
g u is to lowest order constant, independent of u, while £(it) can be replaced 
by 5(s) plus a linear term in s — u. So long as t remains in this linear regime, 
we continue to have a behavior in which the integral grows as log t However, 
as t gets larger, and hence the range of u gets larger, the denominator ceases 
to be a linear function toward the upper limit of the integral. In fact, for 
u — s equal to or of the order of — s, quadratic and higher terms in the 
expansion of £(s + (u — s)) in u — s increase the size of the denominator and 
hence cut off the integral. As that happens the integral increases only very 
slowly with t, and we say that the logarithmic singularity has been cut off. 
In this range of t, we can say that gt(j(s)) has become constant, and takes 
on the value 



9th(s)) = «r(7(*)) " / 2 du/Uu) - g u (j(s))] 



JT 



for A~ 2 < s - T < -s 



(31) 




(32) 
(33) 
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Now apply equation 1311 to a somewhat higher value of — T, i.e. s — T of 
order — s. We can safely replace the g's by the constants fj,, at the cost of a 
small error. The result is 

g t ( 7 (s)) = fx(s) - /* 2 du/Uu) - f,(s)] + 0(1/A) 

for s — T = O(-s) (34) 

Expand to first order in do the integral, and find 

gtM*)) = n(s) + 2\fj- + ^(^V^ + °(VA) (35) 

Equation EHl is the basic result of this section and indeed this paper. It is 
accurate whenever the denominator £(u) — g u (l(s)) always has g neligeble in 
comparison to £ through the region in which u is of order t. This in turn will 
be true whenever —t is very small in comparison to unity. Since equation 1351 
demands that the linearized solution be in its large "x" asymptotic range 
we must have the conditions 

A" 2 «-i + s<l (36) 

Equation OSl represents a small deviation for the critical trajectory. As 
such it should be, and is, identical in form to our earlier result (equation 
I23|) for the deviation. The only change is the the unknown parameter A in 
the earlier result has been replaced by the known parameter fi(s). For all 
values of t near the critical time, — t <C 1, the final term in equation 1351 is 
negligible compared to the first. The second term is large when 

- t » (- S )/V(l-/3) 



5.3 Evaluations of the trace. 



If —to is much smaller than one, equation 1^3 is accurate and one can get the 
value of the trace simply by setting t = to where <7t(7(s)) = t(s). We thus 
find 

j(s) = M (,)[l + Od-to) 1 -^)] + fc^^ + 0(1/ A) 

Here the symbol O indicates real correction terms. Our goal is to find 
the shape of the trace. We use this expression to form (7(5) — 7(0)), using 
fj,(0) = 0, and then take the real and imaginary part of the result, to find 



Re( 7 (s) - 7(0)) ~ £(*) 



JL 

2vr 



Im7(s) 



0/(1-0) 



(37) 
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This result can be extended to values of —to of order unity or larger. 
The mapping that extends the trajectories further toward negative t values is 
conformal and hence differ entiable throughout the upper half plane. As such 
it modifies any short piece of curve, e.g. the near-critical trace, by displacing 
it, magnifying it by the constant factor, C, and rotating it. Because the 
s = point on the trace sits on the real axis, the displacement is parallel to 
the axis and the rotation is nil. Consequently for earlier values of the initial 
time, to , equation 1571 is replaced by 

Re( 7 (s) - 7 (0)) « [/? Im 7 (s)f « Cf (s) (38) 
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Figure 5: This figure plots the real and imaginary parts of the trace against 
the time interval to the singularity. The upper set of points is the real 
part; the lower the imaginary part. For comparison the curves give the 
lowest order asymptotic estimates for the trace from equations I391 The case 
depicted is the same as the one shown in figure ^ There are two adjustable 
constants in the fit. First a constant, with the value 0.345 is subtracted 
from the numerically determined trace, so as to give an estimate of the 
deviation of the trace from its value at zero. Next, the asymptotic curves 
are both multiplied by a constant with value 1.05 so as to make the scale of 
the numerical trace the same as that of the asymptotic one. 
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5.4 Check against numerical results 

To see how well this theory works, we have calculated the orbits, and from 
these the trace for the case in which £(i) = 8 * (1 — tf and to = 0. The 
general shape of the trace was shown in Figure ^ above. 

To check the asymptotic analysis, figure [5] plots the real and imaginary 
parts of the trace against the time to the singularity. The curves on this 
figure are the lowest order estimates of these singularities 

Rer (s) = i(s) + (A(s))- 1 log(A( S ) 2 | S |/2) (39a) 

Imr (s) = 2vriA(s) (39b) 

The data is expected to fit the theory toward the left hand side of this 
figure, and indeed it does so. 

However, log log plots like the one in figureEJcan often conceal difficulties. 
For example the second term in equation I39al is essentially invisible in figure 
|SJ To make a more incisive comparison, take the numerically determined 
7(5) and argue that it can be fit by the asymptotic estimate Tq of equation 
1391 in the form 

7(s) - a = C T (s) + d (T (s)) 2 (40) 

Here a, C and d are real fitting constants. The value of the trace at the 
singular point is given by a. The magnification which might occur during 
the time well before the singularity is represented by the factor C. A further 
non-linear magnification is given by d. A sharp comparison is obtained by 
taking the real and imaginary parts of equation 1321 and diving the left hand 
sides by the right hand sides. The result is shown in figure An agreement 
between theory and numerics would be signaled by a set of points which 
approached unity as the time approached the value at the singularity. The 
curves shown seem to support the asymptotic theory. 
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Figure 6: Compensated traces. This takes the same data as shown in Figure 
El and plots then in such a way as to emphasize the deviations from the 
theory. The real and imaginary parts of the calculated trace (minus the 
value at s = 0) is compared with the theory as described in equation 021 
The lower curve describes Re(7(s) — 7(0)) ; the upper is a plot involving 
Im(7(s)). If theory is right, then the curves should asymptote to unity as 
time goes to zero. It does so in a satisfying fashion. The fitting constants 
are a = 7(0) = 0.345, C = 1.05 and d = 0.05. 



Appendix: Computational Methods 

Two numerical schemes are employed in this paper. One is the standard 
second order Runge-Kutta method [2] which we abbreviate as RK2. The 
other method is a second order accurate version of the well-known [3] square 
root method, abbreviated, which uses 

(g tf (z) - f) 2 = (g u (z) - 2 + 4(t f - U) + 0(t f - U) 3 (41a) 

The parameter, £, is the average forcing, corrected to gain second order 
accuracy 

m f +u)/2) + ^ [§(</) - (41b) 
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To derive equation 141 al start from the exact relation: 

Integrate this between t\ and t2, using the trapezoidal approximation to 
estimate the right hand side. One finds: 



g 2 -2tg-4t\ 2 



71 6 +92 &)(*! - t 2 ) 



One then expands the £'s about the average time to obtain the result in 
equation I41al 

The step size is set by a parameter. STEP, with the actual step size 
being 

(t f -U) = 8TEP\Z-g\/(dZ/dt) (42) 

with some cutoffs to prevent the actual difference, (tf — ti), from getting too 
large or too small. 
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